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Abstract. We provide the first direct lifting of the mass/anisotropy degeneracy for a cluster 
of galaxies, by jointly fitting the line of sight velocity dispersion and kurtosis profiles of the 
Coma cluster, assuming an NFW tracer density profile, a generalized-NFW dark matter profile 
and a constant anisotropy profile. We find that the orbits in Coma must be quasi-isotropic, 
and find a mass consistent with previous analyses, but a concentration parameter 50% higher 
than expected in cosmological A''-body simulations. We then test the accuracy of our method 
on realistic non-spherical systems with substructure and streaming motions, by applying it to 
the ten most massive structures in a cosmological A^-body simulation. We find that our method 
yields fairly accurate results on average (within 20%), although with a wide variation (factor 
1.7 at la) for the concentration parameter, with decreased accuracy and efficiency when the 
projected mean velocity is not constant with radius. 



1. Introduction 

The kincmatical ("Jeans") analyses of near-spherical structures using the spherical 
Jeans equation for stationary systems: 

d , 2^ 2/3 2 GM 
dr r r-^ 

where v is the 3D density distribution of the tracer population, suffer from the fact that 
this one equation involves two unknowns: the mass distribution M(r) and the velocity 
anisotropy (5 — \ — /of. (we make use of the symmetry of spherical systems yielding 
= CTg). In some cases of simple anisotropy profiles /3(r), the Jeans equation can 
be inverted to yield the radial dispersion as a single integral over radii of vM times 
some kernel, and one can insert this solution into the equation that links the line-of- 
sight velocity dispersion to the radial one, to obtain after some algebra, the line-of-sight 
dispersion as a single integral over vAI and some other kernel. 

In what follows, we attempt to lift the mass/anisotropy degeneracy by adding a second 
equation, namely the 4th order Jeans equation for constant anisotropy (Lokas 2002) 

A(^^) + M^^ + 3^^2d£^Q^ 
dr r dr 

where we express the line of sight kurtosis as a double integral over vM and M and a 
third kernel (the details are given in Lokas & Mamon 2003). 



2. Kinematic analysis of the Coma cluster 

We apply our formalism to the nearby Coma cluster, for which we've extracted from 
NED the positions and line-of-sight velocities of 967 galaxies, among which 355 early- type 
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Figure 1. Probability contours (1, 2, and 3 a) for joint fits of crios and (logKios)^ (tlie two 
are found to be independent). Circles indicate the best-fitting parameters. The mass is in units 
of 10^5 Mo . 



galaxies whose kinematics we analyze here (we omit spiral galaxies, which are expected 
to have a non negligible mean radial velocity caused by their infall and subsequent trans- 
formation at cluster pericenter). Our mass model includes 1) a stellar component, which 
we find is well fit by an NFW (Navarro et al. 1996) model; 2) a dark matter component, 
which we assume to be a generalized NFW model with arbitrary inner slope: 

p(r)«(r/a)-"(l + r/a)"-3 , (2.1) 

where a is the scale radius of the dark matter (in general different from that of the 
luminous distribution), and a = 1 for the NFW model. The profiles differ by their inner 
slope —a but have a common outer limiting behavior of r~^; and 3) an isothermal gas 
component, taken from the ROSAT X-ray observations by Briel et al. (1992). 

Using 7 radial bins at i? < 80', we jointly fit our four parameters, mass My within the 
virial radius , velocity anisotropy /3, inner slope a of the dark matter density profile, and 
dark matter concentration c = r^/a. Contrary to the case when only velocity dispersion 
is studied, the minimization procedure now converges (except for a concentration / inner 
slope degeneracy) . The best fits are shown in Figure ^ The best fit parameters are 
My — 1.2 X IO^^Mq (corresponding to dark matter virial radius = 92' = 2.7 Mpc), 
P — —0.13, a = and c = 19 with /N = 6.1/10. In other words, our best fit models 
are isotropic and the generalized dark matter model of equation (|2.1|l provides excellent 
fits to the data. 

Our best fit radial profiles for each half-integer value of a are shown in Figure [21 
The best fit density profiles are indistinguishable beyond 0.03 r^,, and the outer mass 
profiles are very close to the infall estimates of Geller, Diaferio & Kurtz (1999). The 
density profiles are 50% more concentrated than in structures of comparable mass within 
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Figure 2. Best fitting dark matter profiles witfi different inner slopes {left panel) and their 
associated mass-to-light ratios {right panel), for different a increasing upwards from to 1.5 
(with c = 19, 14, 9 and 5, respectively). The dark matter virial radii are r^ = 2.7 Mpc for all 
models shown. The heavy solid curve in the left plot indicates the best NFW model (c = 6) 
from cosmological A^'-body simulations, which is still (marginally) consistent with the data. 

cosmological A^-body simulations in a ACDM Universe (Bullock et al. 2001). Interestingly, 
the M/L profiles have nearly constant slope for the NFW case {a = 1). 

Our analysis differs from previous kinematical analyses (e.g. Merritt 1987; den Hartog 
& Katgert 1996; Carlberg et al. 1997; van der Marel et al. 2000; Biviano & Girardi 2003; 
Katgert et al. 2004), as 1) we have, for a single cluster, a larger sample of galaxies, 
which, given their early morphological type, should be in dynamical equilibrium in the 
cluster potential; 2) we remove pairs from the computation of the velocity moments; 3) 
we include kurtosis in the analysis; 4) we model the dark matter distribution using a 
generalized formula inspired by the results of cosmological iV-body simulations; 5) we 
include hot gas. 

In comparison to studies based upon stacking of many clusters, our analysis of the 
Coma cluster benefits from not having to introduce errors in any stacking procedure, 
and from a cleaner removal of interlopers. On the other hand, the analyses of stacked 
clusters have the advantage of averaging out particular inhomogeneities of individual 
clusters such as Coma, expected in hierarchical scenarios of structure formation and 
observed in Coma, both in projected space (Fitchett & Webster 1987; Mellier et al. 1988; 
Briel et al. 1992) and velocity space (CoUess & Dunn 1996; Biviano et al. 1996). 



3. Kinematic analysis of simulated clusters 

In our Jeans analysis of the Coma cluster, we assumed that the cluster was spherical, 
with no substructure nor streaming motions, and that projection effects were not serious. 
We test (details in Sanchis, Lokas & Mamon 2004) these assumptions using an ACDM 
cosmological A^-body (tree-code) simulation, run by Ninin (1999), with 256'^ particles, 
and with f]„ = 1/3, Qx = 2/3, h = 2/3, and erg = 0.88 (details in Hatton et al. 2003). 

We analyze each of the 10 most massive halos along 3 axes: the major axis, an or- 
thogonal one and an intermediate one. Figure |31 shows the radial profiles of line-of-sight 
velocity moments (mean, dispersion, skewness and kurtosis) along the principal axis. 
These profiles display statistically significant radial variations caused by substructures 
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Figure 3. Projected velocity moments of the dark matter particles in halo 1 measured along 
the principal axis. The upper left panel shows the mean line-of-sight velocity with respect to 

the velocity of the center of the halo in units of the circular velocity at rioo. The upper right 
panel gives the linc-of-siglit velocity dispersion in the same units. The two lower panels give the 
skcwncss {left) and kurtosis (right). In each panel the solid line shows results for particles lying 
inside the sphere of radius rioo, while the dashed line is for all particles. The dotted curve shows 
the fits obtained from the Jeans equations. 
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Figure 4. Fitted (on projected data) values (using 40 particles per bin) of virial mass Mioo, 
concentration parameter c and anisotropy (3 of the ten halos for the three directions of observa- 
tion with respect to the major axis of each halo: 0° (circles), 45° [triangles) and 90° (squares). 
The values measured on 3D data are shown with crosses. 
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Particles AlogA/ioo 


A log c 


Al3 


Alog(crr/cre) 


per bin mean cr 


mean a 


mean a 


mean a 


All -0.03 0.09 


0.20 0.18 


-0.78 1.04 


-0.12 0.12 


40 -0.07 0.10 


0.08 0.24 


-0.20 0.48 


-0.04 0.11 



Table 1. Results of the fitting procedure 



with radial streaming motions. Moreover, the moments vary considerably from halo to 
halo and for the 3 different projections of a given halo. Note also that the projected mean 
velocity differs substantially between all the particles in the projected cylinder and the 
subset within the virial sphere. 

FigureElshows the results of our fits from our Jeans analysis {open symbols) in compari- 
son with our fits to the 3D data (crosses) and TableHsummarizes the statistical accuracy 
of our fitting procedure. The virial mass is typically underestimated by 15 ± 26%, while 
the concentration parameter is overestimated by 20 ± 74% and the ratio Orj <^d under- 
estimated by 9% ± 29%. Hence, our kinematical analysis yields good results on average, 
but with a fairly large dispersion in accuracies (with the most inaccurate estimates for 
halos with non zero mean line of sight velocities at some radii). Halo 4, which departs 
the most from the 3D measures of concentration and anisotropy, happens to have im- 
portant variations in its projected mean velocity profile, whereas nearly all other halos 
have normal mean velocity profiles. This suggests that the joint dispersion/kurtosis Jeans 
analysis is adequate for structures with fairly constant projected mean velocity profiles. 
Given the large dispersion on the concentration parameter, our high value of c for Coma 
is now consistent with the 1.5 times lower value expected from cosmological iV-body 
simulations. 
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